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We investigate imaging of the spatial spin distribution of atoms in optical lattices using non- 
resonant light scattering. We demonstrate how scattering spatially correlated light from the atoms 
can result in spin state images with enhanced spatial resolution. Furthermore, we show how using 
spatially correlated light can lead to direct measurement of the spatial correlations of the atomic 
spin distribution. 
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Experiments with ultracold atom.s in optical lattices 
have become a integral part of investigations into the 
fundamentals of quantum mechanics, precision measure- 
ment and quantum computing, and have become a valu- 
able tool for the simulation of condensed matter prob- 
lems [l|. These experiments typically use time of flight 
absorption imaging 0] , or more recently noise correlation 
imaging Q to examine the atomic state. These methods 
destructively probe the momentum state of the atomic 
cloud with resonant light after releasing the atoms from 
the lattice and allowing them to expand. 

As experiments with optical lattices become more ad- 
vanced, atomic spin has become an important variable 
[J| , and new ways to probe it must be investigated |5n^] . 
Here we consider imaging the spatial atomic spin distri- 
bution non-destructively while the atoms remain in situ. 
In particular we look at ofF-resonant light scattered from 
the atoms, which is then collected by a microscope to 
form an image on a detector. We consider incoming light 
from coherent beams and also from spatially correlated 
beams. Where in the later case we build on the ideas de- 
veloped in the field of quantum imaging [9|, |lO| to show 
how the resolution of images can be enhanced and how 
spatial correlation functions of the atomic spin distribu- 
tion can be measured. These concepts become increasing 
useful for optical lattices where the lattice spacing is often 
similar to the wavelength of the probe light and coherent 
imaging becomes less useful due to the diffraction limit. 

Our approach is organized as follows. We begin by 
describing the interaction between atomic spins and off- 
resonant light that our imaging proposal is based on. 
We then develop a description of how multi-photon light 
states are imaged following their interaction with the 
atoms. We then apply our theory to imaging of a num- 
ber of possible spin states. Finally we show how using 
spatially correlated imaging light can also lead to direct 
measurement of spatial correlation functions. 



I. INTERACTION BETWEEN LIGHT AND 
ATOMIC SPINS 

Many ultracold atom experiments are done with alkali 
atoms and we thus restrict our attention to the inter- 



action between alkali atoms in their ground electronic 
state and detuned light. In the ofF-resonant limit, i.e., 
where the detuning from the 5*1/2 ~^ ^1/21-^3/2 tran- 
sitions is larger than the natural line widths and also 
larger than the Rabi frequency of any particular radiation 
mode, the excited states can be adiabatically eliminated 
from the interaction Hamiltonian. In the dipole and ro- 
tating wave approximations the interaction between the 
radiation and atomic fields is then given by the effective 
Hamiltonian 

m\ ,7712 

^^_. (1) 
where ami,m2 is the polarizability tensor [a,llllll4|- Here 

the atomic field operators ^\^{y) and ^m(r) create and 
destroy atoms in the ground hyperfine state \Fm), and 
E+(r) and E~ (r) are the slowly varying positive and neg- 
ative frequency components of the electric field. 

As discussed in Hammerer et al. 11], the polarizability 
tensor can be decomposed into three irreducible compo- 
nents of rank zero, one and two. The rank zero term 
leads to a spin-independent interaction, which does not 
play a part in the imaging interactions we consider in this 
work. The rank one term is responsible for the Faraday 
effect and is the interaction we use for imaging. The rank 
two term also leads to a spin dependent interaction, but 
is typically at least an order of magnitude smaller than 
the second term and asymptotically goes to zero as the 
detuning increases, and we neglect it here. 

Keeping only the rank one component the interaction 
becomes 

iJ = ai (A) ^ jdvp{v) ■ (E- (r) x E+ (r)) . (2) 



,n,,™. *U(r)(Fm2|F/;i|i^mi)^,„,(r) is 



Here /5(r) = ^ 

a vector operating on the atomic spin and ai(A) is the 

detuning dependent coupling constant (see [11]). 

This effective interaction Hamiltonian describes light 
scattering from the incoming beam, where photons scat- 
ter from the atoms in the lattice and take away infor- 
mation about the spin distribution. These photons can 
then be imaged onto a detector as shown in Figure [T] 



i f, 


f, 


f. 


f. 















FIG. 1. Imaging scheme. Light propagating in the z direction 
interacts with atoms in an optical lattice, then passes through 
a diffraction limited two lens system and is detected at the 
focal plane of the second lens. 



To make calculations for this system practical however, 
the full three dimensional, multimode interaction Hamil- 
tonian must first be simplified. In considering this type 
of imaging our primary concern is the spatial descrip- 
tion of the light as it interacts with the atoms in the 
optical lattice and propagates through the system. With 
this in mind, we assume the photons in the light field 
are approximately monochromatic with frequency w. We 
also consider one-dimensional optical lattices in the x- 
direction, with the result that the useful variation in an 
image is also in the x-direction. The system can then be 
approximately described by a light field that only varies 
in two dimensions, x and z. Numerical calculations of 
images for three dimensional systems show that this sim- 
plified theory captures the nature of the imaging process 
well, and provides a good model for images integrated 
over the j/-direction. 

Tsang [ijll describes the quantization of a two dimen- 
sional electromagnetic field in the monochromatic ap- 
proximation. We generalize Tsang's formalism to include 
both possible polarizations, in which case the electric 
field operator becomes 



E+(x, 
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iKX-]-i^/k''^ —K-^z 



(3) 
where rj = hk/{2eoLyT), j{k) ^ {1 - k^/P)-^/'^, k ^ 
10 1 c = 27r/A, and T and Ly are the normalization time 
scale and y-dimension length scale respectively (see [ij] 
for further details). 
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The operators ai(K) and 62 (k) destroy photons with 
wavevector kx -|- ki/^^K)^ and polarization vectors 
ei(K) = x/7(k)^ — Kz/fc and e2{n) — y respectively. 
They obey the commutator relation [(^^(k), ao-'(K')] — 
S{k — K')5cr,a'- Usiug thcsc operators we can define the 
A'^-photon momentum eigenstate 
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y/n\{N -n)! 
X a\{Ki) . . .a|(K„)4('«n+i) . . .a2(KAr)|0), (4) 

which has n photons with ci-polarization and N — n pho- 
tons with e2-polarization. These states can then be used 
as a basis to express a general state \n)i\N — n)2, which 
again has which has n photons with ei-polarization and 
N — n photons with e2-polarization. 



\n)i\N-n)-^ 



dKl . . . dKN(l){Kl, . . . , Kn) 



-k 



X |ki, . . . K„)i (g) |k„+1, . . . , Km)2, (5) 



where we refer to 4){ni^ . . . , hn) as the momentum distri- 
bution of the A''-photon state. 

We suppose the atoms are in an optical lattice in the 
x-direction with uniform site separation a and expand 
the atomic field operators in terms of lowest band Wan- 

nier functions 'l'm(r) = X^i ^j\ni^(''j)' which we assume 
are not dependent on the spin state of the atom. The 
operator hj^m destroys an atom with spin m at lattice 
site j, centered at r^ — (a;j,0,0). We assume the lattice 
is strong enough so that the Wannier functions have neg- 
ligible overlap between neighboring sites, in which case 
the spin operator becomes 

p(r)= Y. \w{v,)\%^i,,^AFm2W/h\Fm,). (6) 



Integrating over y and z then gives a vector operator 
dependent only on x 



p{x) = / dydzp(v) 



(7) 



In terms of this spin operator and the two-dimensional 
electric field the effective interaction Hamiltonian be- 



comes 



H^'^li^ldx 
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In deriving Equation ([S]) we used the approximation which is very good for the imaging scenarios we con- 
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2 l{y/k^^-^/W^)z 



dydz\w{rj 
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sider, where the Wannier functions are confined in the 
z-direction tightly compared to the variation of the elec- 
tric field in the z-direction. 

Furthermore, for the imaging parameters we consider, 
the term proportional to p{x) ■ z — Pz{x) in the effec- 
tive Hamiltonian is much larger than the others and de- 
termines the dominant features of the images. In the 
following we neglect the other terms in our theoretical 
description. These terms lead to small changes that can 
be easily calculated and are included numerically in the 
calculations used to produce our figures. 



II. LIGHT SCATTERING AND N-PHOTON 
IMAGING 

Having derived a description of the interaction between 
the atomic spins and the light field, we now consider 
how this interaction leads to light scattering and how 
the scattered light can be used to image the spin dis- 
tribution. Taking the atomic spin quantization axis to 
be z, the spin operator in the 2-direction is diagonal in 
atomic spin, that is Pz{x) = '^j Pz,j\w{x — Xj)^, where 

Pz,j = Sm'^^j,™ and fij^rn = b]^mP],m- Thc atomic evo- 
lution in the optical lattice due to intersite tunneling and 
atom-atom interactions occurs on a time scale of order 
the inverse recoil frequency [IJ . We will assume that the 
duration of the light matter interaction is short compared 
to this time scale and neglect the atomic evolution dur- 
ing the interaction. First order perturbation then mixes 
the original state of the atom-photon system |-0) with 



dx / dKidK2e*('"^-'"i)^/5^(a;) 



X f ^al(Ati)a2('*2) - ^-p\al{n^)d,{K^)) \4>) (10) 
V7('«i) 7(k2) / 

This first order scattering process will trans- 
form the initial momentum distribution of the pho- 
tons (f>i{Ki, . . . , kn) into the scattered distribution 
(J)s{ki, . . . , Kjv). The imaging system in Figure[T]will then 
map the scattered momentum distribution to a new mo- 
mentum distribution at the detector ^^(ki, . . . , km)- Due 
to the finite extent of the lenses in the imaging system, 
the transverse momentum of the photons reaching the de- 
tector will be restricted to a finite bandwidth, such that 
|fi:| < K; = fcsin^ where sm.9 is the numerical aperture of 
the optical system J15l.[l6|. This limit, which corresponds 
to the Rayleigh-Abbe diffraction limit, requires 



4>d(,Ki, . . . ,kn) 



if any |k„| > ki 

0s (ki, . . . , kn) otherwise 

(11) 
where for simplicity we have let /i = /2- 

We now have a description of the A^-photon state at the 
detector and we are left with the problem of extracting 
useful information from the state. A measurement of 
the intensity of the light at the detector provides the 



simplest form of image, but does not make full use of 
the multiphoton state. Using an A'^-photon absorbing 
material that creates a signal when all A^ photons arrive 
at the same point allows us to make better use of the 
iV-p hoton nature, but may be subjectjx) low efficiency 
[13| . Recently it was demonstrated in [15] that intensity 
measurements at the detector followed by post-processing 
can reveal similar information to A^-photon absorption. 
As this method does not rely on the photons all being 
coincident at a single point it has much higher efficiency. 
The variable required in these measurements is 
the A^-photon intensity distribution, I{xi, . . . ,xn) ~ 

J:a,,...,ajE-JxN) ■ . .-B- (Xi)^+ (Xi) . ..E+JXJ,)). 

At the detector for the state \n)x\N — n)y this is given 
by 



Hxi,. 



,Xn) 



^77^|Vd(P(xi,...,x^))|' (12) 



where the sum is over all coordinate permutations P and 



ipd{xi,. . . , xn) = dni... duNlini) ■ . . i{kn) 

1 f ^ \ 

is a A^-dimensional transform of the momentum dis- 
tribution. The transform ipi{xi, . . . ,xn) of the input 
light is defined in the same way. Note that if all pho- 
tons are in the same polarization state, the momen- 
tum distribution must obey bosonic symmetry so that 
(/)(..., K„, ... , Km, ■■■) = (/>{■■■, Km, . . . , K„, . . .) for ah n 
and TO. In which case ip{xi, . . . , xn) must obey the same 
symmetry and is invariant under coordinate permuta- 
tions. 



III. COHERENT SPIN IMAGING 

Before studying multiphoton imaging, we first look at 
the image created due to coherent laser photons scat- 
tering off the atoms in the optical lattice. We suppose 
the photons enter the system with e2-polarization and 
momentum distribution 4>i{K). First order scattering ac- 
cording to Equation ()10p will then result in ci-polarized 
photons with a momentum distribution that depends on 
the atomic field operators according to 



(^s(k) 
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7(k) 



dn' / dx7(K')e'^'''""^"0»(K')/5.(a;)- 



(14) 

Here, because light scattering off the atoms results in an 
orthogonal polarization, the unscattered light can be fil- 
tered out and we can detect an image due only to the 
scattered light [6]. The intensity distribution at the de- 
tector will then be 



I{x) 



ri{ip\{x)'il)d{x)) 



(15) 



dK,'y{K)(j)s{K)e^ 



2k; dx'sinc{Ki{x — x'))ipi{x')pz{x'). (16) 



where 



The result is a diffraction hmited image of the spin den- 
sity weighted by the photons original spatial distribution. 
The integration over the sine function in the image ampli- 
tude ipdix) blurs the image, but has the useful side effect 
that the image depends on the spin correlation function 
{p.{x')p.ix")). 



IV. TWO-PHOTON SPIN IMAGE 

We next consider the scattering of a two-photon state 
from the atoms in the optical lattice and the result- 
ing two-photon intensity distribution. We take the ini- 
tial state to have one ei-polarized photon and one €2- 
polarized photon, then first order scattering from this 
state will result in a state with two photons with the 
same polarization. The photons in the final state will 
then be indistinguishable and knowledge of which photon 
was scattered is lost. Using polarization filters or beam 
splitters, we can then detect only two-photon states with 
two photons in the same polarization state. The result- 
ing image is then due to pairs of photons containing one 
scattered photon and one unscattered photon. If we fil- 
ter out €2 -polarized photons, the two-photon intensity 
distribution is 

I{XUX2) - {E+{xi)E+{x2)E^{x2)E^{xi)) 

^2rf{ip\{xi,X2)ipd{xi,X2)) (17) 

where 

V..(x.,..)^V2.,/d.AW(.inc(.,(.,-..))*(x..) 

+smc{Ki{x2 - x))'il!,{xi,x)) . (18) 

The amplitude 'ipdixi,X2) consists of two terms due to 
the two indistinguishable paths the photons could have 
taken to reach xi and X2, that is, the scattered photon 
reaching xi and the unscattered photon reaching X2 and 
vice versa. For simplicity here we have assumed that the 
photon source is also limited in transverse momentum by 

Kl- 

The use of two-photon absorption imaging can now 
lead to higher resolution than that of the coherent spin 
image discussed in the previous section. If the initial two- 
photon state is generated by parametric down conversion, 
then the transverse momenta of the photons is typically 
anti-correlated [17|. If we take the initial momentum 
distribution to be approximately (piini, K2) = 6{ki + K2) 
for |ki|,|k2| < Ki and zero elsewhere, then the spatial 



distribution is ipi{x,x') = —smc{Ki{x — x')). In this case 
a two-photon absorption image depends on the amplitude 



ipdixi,xi) 



VSk 



dxpz{x)amc^{Ki{xi — xj). (19) 



Comparing this with the coherent imaging amplitude in 
Equation (J16p . we see that using two photons here re- 
sults in the sine function being raised to the power of 
two, allowing an image with higher resolution. In gen- 
eral for an iV-photon state this type of imaging results in 
an improvement in resolution that scales like ^/N [9|. 

A two-photon absorption image requires two photon 
detectors. A more efficient way to acquire the resolution 
improvement is to calculate the centroid distribution of 
the two-photon intensity distribution [1^ . In place of the 
Cartesian coordinates we can express the system in terms 
of the centroid and relative positions of the two photons 



X 



Xi +X2 



e 



Xi - X2 



(20) 



For the parametric down converted source described 
above the centroid probability distribution (unnormal- 
ized) is 



d^i^Ux + ^,x- OMx + LX- 0) 

8(«/)4 



dxdx'{p^{x)p,_{x'))K{X, X, x'), (21) 



where 



K{X,x,x')^ / d(smc{Ki{X+^-x))smc{Ki{X-^-x)) 

X smc{Ki{X + C - x'))smc{Ki{X - ^ - x')). (22) 

This results in a similar resolution improvement as for 
the two-photon absorption above, but is expected to be 
much simpler to implement in an experiment. 

If we could measure the amplitude (ipixi , X2)) directly, 
for example by spatial quantum state tomography jl8| . 
then the centroid measurement applied to the amplitude 
would result in an image with higher resolution still. Us- 
ing the identity 



dxsmc{b{x + y))smc{b{z — x)) = —smc{b{y + z)) 

(23) 
the probability distribution (unnormalized) of the cen- 
troid becomes 



d^mx + tx-0) 



dx{pz{x))smc(2Ki(X — x)) 



(24) 



The factor of two is now inside the sine function, lead- 
ing to a factor of two resolution improvement over the 






FIG. 2. Normalized images of a lattice containing atoms all 
in the m = 1 state with a defect, a) The lattice spin density 
to be imaged, b) The coherent image, c) The two-photon 
absorption image, d) The centroid distribution for the two- 
photon intensity distribution (Eq. ((21}). e) The centroid of 
the two-photon amplitude (Eq. (I24p ). Here the lattice site 
separation is 0.9A and sin6' = 2/3. 



coherent case. However, the requirement for quantum 
state tomography means that this would be significantly 
more difficult to achieve in an experiment. 

In Figure [2] we compare the various possible images de- 
scribed above. The images are of an optical lattice with 
a single m — 1 atom at each site except for the missing or 
TO = atom at a; = — 0.45A. The lattice site separation 
is 0.9 times the wavelength of the imaging light and the 
optical system has a numerical aperture of 2/3. For these 
parameters the coherent image allows the identification 
of the defect in the optical lattice filling, but does not 
resolve the individual sites. The two-photon absorption 
image and the centroid of the two-photon intensity dis- 
tribution both resolve the individual sites. The centroid 
of the two-photon amplitude resolves the sites with an 
improved visibility. 

The precise resolving power of each setup depends on 
the lattice properties, but a rough guide it that in the 
coherent image it becomes difficult to resolve individual 
sites for site spacings below a ~ A/ sin 6*, while for the 
two-photon amplitude centroid this limit is halved. For 
the two-photon absorption and two-photon intensity cen- 
troid images individual peaks associated with each lattice 
also develop around a ~ A/2 sin as in the amplitude 
centroid image, but the visibility of these peaks is lower 



than for the amplitude centroid image. 



V. IMAGES OF EXAMPLE SPIN STATES 

We now apply our imaging proposal to example spin 
states of atoms in an optical lattice. For an optical lattice 
with a single spin-1 atom at each site, the ground state 
potentially consists of dimers [l^- This occurs when 
neighboring sites form a spin singlet (total spin equal 
to zero). Denoting site j's spin states m — +1,0,-1 by 
|+)j,|0)j and \—)j respectively, a dimerized state of an 
A'f-site lattice (M even) is represented by 

Af/2 

-1-10)2,10)2,-1) 
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(^^(|+)2,|-)2,-l + |->2,|+)2,- 

(25) 
The coherent imaging intensity can be rewritten in 
terms of the site correlation functions as 

I{x) = 4f]Kf y^(pz,jPz,r)f{x - Xj)f{x - Xr) (26) 



j,r 



where 



f{x)— / dx'\w{x')\'^smc{Ki{x — x')), 



(27) 



For the dimer case 



and |w(a;)p = J dydz\w{x,y, z)\'^ 

{Pz,jPz,r) = otherwise. The intensity then becomes 



{pl.j) = 2/3, {pz,2j-iPz,2j) = -2/3 for integer j and 



I{x) 



8r]K 



2 M/2 

E(/(^ 



X2j-l)^f{x^X2j))\ (28) 



From this form we see that the atomic correlations lead 
to interference in the resulting image. In particular, at 
the centre of each dimer pair, x = {x2j-i + X2j)/2, the 
contribution to the image from the atoms in the dimmer 
at X2J-1 and X2j cancel, and we are left with a back- 
ground intensity due to the other atoms. This cancel- 
lation leads to local minima in the image at the cen- 
ter of each dimer pair. The cancellation does not occur 
at a; = {x2j + X2j+i)/2 between atoms from separate 
dimers. For a lattice spacing a < 1/(2 sin 6*) the images 
of sites 2j and 2j -I- 1 are no longer resolved and merge 
at a; = {x2j + 0:2,-1-1 )/2 to form local maxima. Thus be- 
low the diffraction limit, the image is an oscillation with 
double the spatial period of the lattice due to the par- 
ing of atoms. This pattern persists until oscillations with 
double the lattice period can no longer be resolved. 

Another potential state of the spin-1 chain is a trimer- 
ized state, where three neighboring sites form a singlet 
|20| . A trimerized state of an Af -site lattice (M multiple 
of three) is 



J\//3 



1 



QS)-^(l + o-), + |-+o), + |o-+), 

-| + -0),-|0 + -),-|-0-H),) (29) 



where \m1m2m2) 



|mi)3j_2|m2)3j-i|m3)3j. In this 



case (pIj) = 2/3, {Pz,3j-2Pz,3j-l) == {Pz,3j-2Pz,3]) 

{pz,3j-iPz,3j) = -2/6 for integer j and {pz,jPz,r) = 
otherwise. 

The coherent image intensity then becomes 



I{x) 



irjnf 



M/3 



o XI ((•^(^^^3:'-2)-/(a;-a;3j_i))2 

+ (/(a; - 2;3j-2) - f{x - xsj))"^ 

+ {f{x~X3j^,)^f{x-X3,)f). (30) 



As in the case of the dimer, the atomic correlations lead 
to interference in the image. In this case the interfer- 
ence results in an oscillation with three times the spatial 
period of the lattice, with local minima at the center of 
each trimer and maxima at the points between separate 
trimers. 

The same effects occur in the two-photon images of 
both the dimer and trimer, but due to the greater resolv- 
ing power they are visible for smaller lattice site spac- 
ings. In Figure [3] the coherent and two-photon intensity 
centroid images of the the dimer and trimer states are 



compared with the unpolarized state that has {pz 



0, 



(p2 j) = 2/3 and {pzjPz,r) = for j ^ r. The im- 
ages are generated for parameters where the dimer os- 
cillations become difficult to resolve in the coherent im- 
age. The dimer oscillations in the coherent image become 
unresolvable for lattice spacings below a ~ 0.2A/sin6', 
whereas the oscillations in the two photon image become 
unresolvable at a ~ 0.14A/sin6'. The oscillations in the 
trimer state become unresolved in the coherent image be- 
low a '^ 0.16A/ sin0, and in the two-photon image below 
a -0.1 A/ sin 6*. 
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FIG. 3. Normalized coherent (left column) and two-photon 
intensity centroid (right column) images of an 18 atom spin 
chain in the (a-b) unpolarized state , (c-d) dimer state and 
(e-f) trimer state. The images are for a lattice spacing 0.4A 
and sinO — 1/2. 



J dx\'w{x)\'^smc {ki{x ~ a)) then the main contribution 
to the two-photon distribution at (a;i,X2) = {xj,Xj+i) is 

{^pl{xj,Xj+i)'il;d{xj,Xj+i)) « -^{{pz,j +Pz,j+i)'^) 



VI. SPATIAL CORRELATION 
MEASUREMENTS 

Another potential use of two photon imaging is to ex- 
amine local spatial correlations in the atomic lattice that 
cannot be accessed using methods that probe the entire 
sample [3, Q . If we have light that is spatial correlated as 
(j)i{x,x') = -^{smc{Ki{x~x' — a)) + smc{Ki{x — x' + a))), 
so that the two photons are separated by a, then the two- 
photon image amplitude (Equation (|18p) becomes 



■4>d{xi,X2) ^ — dx{pz{x) + pz{x + a)) 



X (sinc(Ki(a;i — a;))sinc(K;(2;2 — x — a)) + 

sinc(K;(a;i — x — a))sinc(K/(x2 — x)). (31) 

A measurement of the two-photon intensity distribu- 
tion at xi = X and X2 = x + a then depends on 
{{pz{x) + pz{x + a)y). For an optical lattice, if the two 
photon resolution is such that J dx\w{x)\'^smc (kix) 3> 



X / dx\w{x)\^sinc^{Kix) . (32) 



This allows us to access the neighbor correlation us- 
ing {{Pz,j + Pz,j+if) = {pI, + /5l,-+i) + 2{pz,jPz,j+i) 



For the dimer state the atoms of a dimer pair have a 
neighbor correlation leading to I{x2j-i,X2j) ~ (pij-i + 
P2j) + 2{p2j-iP2j) = 0. There is no neighbor cor- 
relation between atoms of separate dimers leading to 
I{x2j,X2j+i) « {p%2j + pI,2j+i) = 4/3. For the trimer 
state the correlation between neighbors in a trimer leads 

to I{x3j^2,X3.j-i) = I{x3j^2,X3j) = I{x3j^i,X3j) W 2/3. 

While atoms of separate trimers have no correlation lead- 
ing to /(a:;3-,, a;3-,+i) « 4/3 

In the first row of Figure [4] we compare the resulting 
two-photon images for the trimer and dimer states with 
the unpolarized state. The correlation between atoms 
in the dimer state leads to the removal of every sec- 
ond peak compared to the unpolarized state distribution, 
while the trimer correlations lead to partial cancellation 
of the peaks corresponding to the correlated atoms in 
each trimer. 







FIG. 4. Images of two-photon intensity for spatial correlation measurement using photons initially correlated to be separated 
by a (top row) and 2a (bottom row). Left column: unpolarized state. Middle column: dimer state. Right column: trimer 
state. The images are for sinO — 1/2 and a = X. 



Similarly, we could also consider light correlated so 
that photons in the initial pair were separated by 2a. 
The two-photon image would then probe the correlation 
{Pz,jPz,j+2)- For the dimer no such correlations exist and 
cancellation does not occur as it did in the previous ex- 
ample. The trimer does have correlations over this range 
and partial cancellation of peaks occurs. This behavior 
is shown in the second row of images in Figure |4l 



VII. CONCLUSION 

As experiments with ultracold atoms in optical lat- 
tices become more complex, there is a need for methods 



of probing the resulting states that provide more detailed 
information. Here we have demonstrated how coherent 
in situ imaging can distinguish between spin states in 
an optical lattice. Further we have shown how the use 
of spatially correlated light can improve the resolution 
of these images and allow us to distinguish spin states 
beyond the limits imposed by Rayleigh diffraction. We 
have also demonstrated a method that uses spatially cor- 
related light to directly probe the local spatial correlation 
function of the atoms. While we have restricted most of 
our attention here to two-photon imaging due to the rela- 
tive ease of generating two-photon states in experiments, 
these methods can be readily generalized to higher pho- 
ton number states. 
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